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Abstract. In this paper we prove Garvan's conjectured formula for the square of the modular dis- 

f~~>- [ criminant A as a 3 by 3 Hankel determinant of classical Eisenstein series E2n ■ We then obtain similar 

Cn ' formulas involving minors of Hankel determinants for SarA™, for m = 1,2,3 and r = 2, 3, 4, 5, 7, and 

E14A . We next use Mathematica to discover, and then the standard structure theory of the ring of 

modular forms, to derive the general form of our infinite family of formulas extending the classical for- 

| - ■ mula for A and Garvan's formula for A^ . This general formula expresses the nXn Hankel determinant 

^— ( ' det(i52(i+j)('z))i<i,j<jT, as the product of A"'~^(t), a homogeneous polynomial in _E| and Eq, and if 

(-H , needed, E4,. We also include a simple verification proof of the classical 2 by 2 Hankel determinant 

formula for A. This proof depends upon polynomial properties of elliptic function parameters from 

Jacobi's Fundamenta Nova. The modular forms approach provides a convenient explanation for the 

determinant identities in this paper. 
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O . 1- Introduction 

m 

In this paper we prove Garvan's conjectured formula [11] for the square of the modular discrim- 

^^ . inant A as a 3 by 3 Hankel determinant of classical Eisenstein series i?2n- We then obtain similar 

(<— ^ \ formulas involving minors of Hankel determinants for i?2r A'", for m = 1,2,3 and r = 2, 3, 4, 5, 7, 

C^ • and i?i4A'^. We next use Mathematica [37] to discover, and then the modular forms approach of 

f-| ■ [31, pp. 88-93], as outlined in [5], to derive the general form of our infinite family of formulas 

extending the classical formula for A and Garvan's formula for A^. This general formula expresses 

the n X n Hankel determinant det(£'2(i+j)(Q'))i<i,j<n as the product of A"~^(r), a homogeneous 

polynomial in E^ and Eq, and if needed, E4. We also include a simple verification proof of the 

classical formula for A in (1.5) below. This proof depends upon polynomial properties of elliptic 

function parameters from Jacobi's Fundamenta Nova [16]. The modular forms approach provides 

a convenient explanation for the determinant identities in this paper. 

The modular discriminant A is defined in [2, Entry 12, pp. 326] and [28, Eqn. (6.1.11), pp. 
196] by means of the following definition. 

Definition 1.1. Let q := exp(27rir), where r is in the upper half-plane TC. We then have 

A{r)^A{q):=qll{l-q^r- (1-1) 

r=l 
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The Fourier expansions of the classical Eisenstein series En{T) as given by [2, pp. 318] and [28, 
pp. 194-195] are determined by the following definition. 

Definition 1.2. Let q := exp(27rir), where r is in the upper half-plane TC, and take y := Ini(r) > 0. 
Let n = 1, 2, 3, • • • . We then have 

rq"^ 3 



E,{r) ^ E,{q) : = 1 - 24 J] -^ - ^, (1.2) 

■^ 1 — vry 



^=1 - Q ^y 



and for n > 2, 



E2n{T) ^ E2n{q) ^ = 1 " ^ ^ ^^^^ , (1-3) 



>2n -, 

r=l 



with the i?2n the Bernoulli numbers defined in [7, pp. 48-49] by 

t °° t^ 

^— :=Vi?„-, for|i|<2^. (1.4) 

e* — 1 ■^^ n! 

n=0 

The fundamental classical formula for the modular discriminant A is provided by the following 
theorem. 

Theorem 1.3. Let q := exp(27rir), where r is in the upper half-plane TC. Let A(r) and i?2n = 
E2n[<i) be determined by Definitions 1.1 and 1.2, respectively. Then, for \q\ < 1, 

AiT) = j^{El-Ei). (1.5) 

Early elliptic function references for (1.5) are [14], [15, pp. 561], [18, Eqns. (1) and (2), pp. 
154], [25], and [26, pp. 27]. (Both Hurwitz and Molin replace q by q'^.) All of these authors 
contributions refer to earlier background developments in [8, 9, 29]. The chapter notes in [6, pp. 
95, pp. 136-137] provide an excellent summary of Dedekind's fundamental [8, Eqn. (3), pp. 281; 
Eqn. (13), pp. 283; Eqn. (24), pp. 285], Dedekind's comments in [9] on Riemann's work in [29], 
and the remarks of Fricke [10] and Molin [26, pp. 28] on related methods of Jacobi and Hermite. 
The classical elliptic function methods for proving (1.5) are discussed in [6, pp. 58-72], [13, pp. 
409, pp. 481], and [22, pp. 125-140, pp. 177]. An elementary proof of (1.5) was later found by 
Ramanujan in [27, Sections 5, 7 and 10]. This and additional work of Ramanujan involving (1.5) 
is surveyed in [3, pp. 114-140], [4, pp. 43-50], and [34, pp. 1-18]. Two additional elementary 
proofs of (1.5) are described in [31, pp. 95-96]. Recent references for (1.5) are [2, Entry 12(i), pp. 
326], [6, Theorem 7, pp. 71], [28, Eqn. (6.1.14), pp. 197], [30, Eqn. (8), pp. 55; Theorem 8, pp. 
70], and [31, Eqn (42), pp. 95]. (The work [19] is a very useful introduction to [30].) Additional 
applications of (1.5) also appear in [1, 28, 31]. 

After seeing [23, Theorems 2.1 and 2.2, pp. 15006] and an early version of [24, Theorems 1.5, 
1.6, 5.3-5.6] Garvan [11] observed via the well-known relation E^ = Eg that (1.5) immediately 
becomes the Hankel determinant formula 



Mr) = T7W(^4i^8 - El) = j^ det 



E4 Eq 
Eq Es 



(1.6) 



He then conjectured the following theorem. 
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Theorem 1.4 (Garvan). Let q := exp(27rzr), where r is in the upper half-plane 7i. Let A(r) 
and E2n = I^2n{(l) be determined by Definitions 1.1 and 1.2, respectively. Then, for \q\ < 1, 



A\r) 



691 

(1728)^-250 



det 



Ea Ep. 



Er Ea 



Ea 



E 



10 



Ea E 



10 



E 



12 



(1.7) 



Proof. Substitute the following three weh-known relations from [27, Table I., pp. 141], [28, pp. 
195] into the 3 by 3 Hankel determinant in (1.7). 



Ea — E' 



'4) 



E- 



10 



EaEr 



-C/12 — ■HTTT-C'/i + 691-^6- 



691 



{U 



Simplifying, factoring, and applying (1.5) then gives the A^(r) on the left-hand-side of (1.7). □ 

For A"(r) with n > 2, formulas analogous to (1.6) and (1.7) generally require a suitable n + 1 
by n + 1 determinant on the right-hand-side and an additional polynomial factor in Ef and Eq on 
the left-hand-side. This extra polynomial factor can often be simplified by relations such as (1.8). 

We organize the rest of our paper as follows. In Section 2 we first apply recursive methods 
to obtain our 19 determinental formulas expressing small powers of A(r), multiplied by a single 
Eisenstein series, as a suitable constant times a certain minor of a Hankel determinant of the E2r- 
These formulas were motivated by Ramanujan's consideration of -B2r A, for r = 2, 3, 4, 5, 7, in [27, 
Section 16], and the discussion of these i?2rA in [33, pp. 302]. The minors here were initially 
motivated by the re by n minors of the n -|- 1 by re -|- 1 Hankel determinants as discussed in [17, 
pp. 244-250]. We next use Mathematica [37] to discover, and then the modular forms approach 
of [31, pp. 88-93], as outlined in [5], to derive the general form of our infinite family of formulas 
extending (1.6) and (1.7) that involve A"~^(r) and an re by re Hankel determinant of the E2r- 

In Section 3 we follow Jacobi's analysis in [16, Section 42] and utilize the Fourier series for the 
Jacobi elliptic function ns^ to write down a formula for the Eisenstein series E2n, for re > 2. We 
then apply [16, Eqn. (2.), Section 36] to put together a simple verification proof of the classical 
formula for A in (1.5). 

Symmetry properties of the coefficients in the Maclaurin series expansion of ns^ strongly suggest 
that formulas such as (3.8) and (3.9) in Theorem 3.1 will be useful in a further study of the 
determinental formulas in Section 2. 

2. Additional determinental formulas involving powers of A 

Our 19 determinental formulas in Theorem 2.3 involving small powers of A(r), and the infinite 
families of identities in Theorem 2.5 are partly motivated by the determinants in the following 
definition. 



Definition 2.1. Let {cj^}^^ be a sequence in C^, and let m,n = 1,2,3, 
Xn to be the determinants of re x n square matrices 

i/«^i^(i)(K}):=det 



Xi"^^ = X^H{c.}) := det 



/Cl 


C2 


Cn-1 


Cn \ 




C2 


C3 


Cn 


Cn+1 


; 


\Cn 


Cn+1 • • 


C2n-2 


C2n-1 / 




/Cl 


C2 


C-n — m 


Cn-77i+2 


C2 


C3 


Cn — 771+1 ^n — m.-\-3 


\Cn 


Cn+1 • • 


C2n — 7n — 


1 C2n- 


-m+1 



We take Hn and 



Cn+l\ 

Cn+2 

C2n 



(2.1) 



(2.2) 
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The matrix for Xn is obtained from the matrix for H^j^i by deleting the (n — m + l)-st column 



and the last row. Others denote Xn ^Y Xn- We also have Hq 



1 and X 



(m) 



0, if n < m. 



Formally, x 



(0) 



H, 



(1) 



Applications of the Hankel determinants Hn and determinants Xn to continued fractions and 
orthogonal polynomials are discussed in [17, pp. 244-250]. An excellent survey of the literature 
on Hankel determinants can be found in Krattenthaler's summary in [20, pp. 20-23 ; pp. 46-48]. 

Our eventual aim is to study the combinatorial and geometrical implications of the determinants 

Hn {{E2{u+i){<l)}), Xn {{E2{u+i){(l)]), and additional minors, where -62(^+1) (g), with v>l, are 
the Eisenstein series in Definition 1.2. A starting point is Theorems 2.3 and 2.5 below. 

Motivated by our proof of Theorem 1.4 we first consider Ramanujan's [2, Entry 14, pp. 332] 
recursion for the E2r given by the following theorem. 

Theorem 2.2 (Ramanujan). Let q := exp(27rir), where t is in the upper half-plane TC. Let 
E2n{'T) be determined by Definition 1.2. For convenience, with n > 1, define S2n by 



S- 



2n 



-iB 



2n 



An 



E2n{T). 



Lf n is an even integer exceeding 4 then 



(2.3) 



(n+2)(n+3) q 
2n{n-l) ^^"+2 



20 



n 



SaD. 



4Jn-2 



[(n-2)/4] 



2r 



(n + 3 — 5r) (n — 8 — 5r) 



5ir-2){r + 3)}S2r+2Sn-2 



(2.4) 



where the prime on the summation sign indicates that if (n — 2)/4 is an integer, then the last term 
of the sum is to be multiplied by ^ . 

Substituting n = 6,8, 10 into (2.4) yields the relations in (1.8), and setting n = 12 in (2.4) leads 
to the well-known relation 

Eu = EJEq. (2.5) 

We next utilize (1.5), Theorem 2.2, and Mathematica [37] to derive our determinental formulas 
for £^2r'^™') for m = 1,2,3 and r = 2,3,4,5,7, and E^A'^. In each of the 19 identities below we 
first used (2.4) to write all the i?2r in the determinants on the right-hand-sides as polynomials in 
E4 and Eq. Simplifying, factoring, applying (1.5), and then referring to (1.8) and (2.5) as needed 
yielded the left-hand-side of each identity. We have the following theorem. 

Theorem 2.3. Let q := exp(27rir), where t is in the upper half-plane TC. Let A(t) and i?2n = 
E2n{(l) be determined by Definitions 1.1 and 1.2, respectively. Then, for \q\ < 1, 



^4A(r) 



691 



1728-250 



det 



Ea Er 



hia hj 



12 



(2.6) 



^6A(r) 



691 



1728-250 



det 



Ea E 



12 



h/R hi 



14 



(2.7) 
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E4 Eio 



E^A(t) 



3617 J . 

1728-3-7^-ll ^^^ 



Elo EiQ 



691 



1728-441 



det 



Es EiQ 
Eio E12 



^inA(r) 



691-43867 
1728-2-32-5'*-72-13 



det 



E4 E12 
Eio Eis 



691 



1728-250 



det 



Es E12 
Eio Ei^ 



EmMt) 



691-593-131 
1728-2-3-53-72-11-13 



det 



Ea E- 



691-3617 
1728-2-3-53-72-13 



det 



14 
E12 E22 

E\o Ei4 
E\2 EiQ 



E^A^t) 



69V 

(1728)2-(21)2-250 



det 



E,A^(t) 



691^ 

(1728)2-(250)2 



det 



E,A^(t) 



(691)"'-3617 



(1728)2-23-3-53-72-467 



EioA\t) 



(691)2-3617 
(1728)2-22-3-56-72-13 



E4, Eg Eio 
Eq Eio E12 
Eg E12 Ei4 
Eq Eg, Eio 
Eg, Eio E12 
Eio E12 Ei4 

Eq Es E12 

■ det Eg Eio Ei4 

Eio E12 EiQ 

Eq Eio E12 

det Eg E12 Ei4 

Eio Ei4 EiQ 



(2i 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



i^MA^fr) 



(691)^-3617-43867 
(1728)^ ■2e-3-5^ -7^ ■97-7213 



det 



E^A'(t) 



(691)^-3617 



(1728)3-24-3-56-72-467 



det 



^fiA^fr) 



(691)^-43867 
(1728)3-2^-32-5«-7^-131 



det 



^sA^fr) 



(691)^-(3617)^ 
(1728)3 -24-32 -513 -7-*- 13-467 



E.oAHt) 



(691)-'-3617-43867 
(1728)3-27-3-5'3-72-97-7213 



E,4A'(t) 



(691)3 -(3617)2 -43867-283-617 
(1728)4-23-32-513-72-31-3503110621 *^^ 
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Es Eio Ei4 
Eio E12 EiQ 

E12 £^14 -E'lS 

E4 Eq Eg EiQ 

Eq Eg EiQ E12 

Es Eio E12 Ei4 

Eio E12 Ei4 EiQ 

E4 Eq Eg E12 

Eq Es Eio E14 

Es Eio E12 EiQ 

Eio E12 Ei4 Eis 

E4 Eq Eio E12 

Eq Es E12 Ei4 

Es Eio Ei4 EiQ 

Eio E12 EiQ Eis 

E4 Es Eio E12 

Eq Eio E12 Ei4 

Es E12 Ei4 EiQ 

Eio Ei4 EiQ Eis 

Eq Es Eio E14 

Es Eio E12 EiQ 

Eio E12 Ei4 Eis 

E12 Ei4 EiQ E20 



det 



det 



(2.18) 



(2.19) 



(2.20) 



(2.211 



(2.22) 



(2.23) 
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E4 Eq Eg EiQ Ei4 



EuA'^ir) 



(eOl)"*- (3617)^ ■43867- 131-283-593-617 . , 
(1728)6. 59 •75- 11- 13-67- 257-43721 ^^^ 



Eq Eg, EiQ E12 EiQ 

Eg, EiQ E12 Ei4 Eis 
Ew E12 Ei4 EiQ E20 



E^o El 4 EiF, E- 



^14 



16 



^18 



E. 



22 



(2.24) 



The determinants in (1.6), (1.7), (2.14), (2.15), (2.19)-(2.22), (2.24) are of the form F^^\ H^'^\ 
xf\ xi\ Hl^ , X4, xf , xf, and X5, respectively, with entries c„ = E2{^+i){q). The rest are 
certain other minors of Hn \{E2[v+i){(l)})-, for suitable n. 

The determinant evaluations in (2.6)-(2.24) can also be proven, as pointed out in [5], by the 
methods in the modular forms proof of Theorem 2.5 below. Being able to apply the first part of 
this proof leads to a characterization of a large class of determinants of Eisenstein series which 
have evaluations analogous to those in (2.6)-(2.24) and Theorem 2.5. It turns out we only have to 
consider minors of Hn ({£'2(!/+i)((i')})) for suitable n. We have the following proposition. 

Proposition 2.4. Let A he any n x n square matrix whose entries are Eisenstein series E2r, 
with r > 2. Suppose there are no repeated rows or columns. Recall that the weight of a product 
-E'2ri-E'2r2 ' ' ' -E'2r„ of Eiscnstcin series is 2(ri + r2 + • • • + r„). Then, each term in the n x n 
determinant detA has the same weight if and only if det A is ±1 times some n x n minor of a 
Hankel determinant Hm {{E2[i,+i){q)}) = det(£'2(i+j)('?))i<i,j<m of Eisenstein series, withm > n. 

Proof. First, assume that detvi is some n x n minor of det(£'2(i+j)('?))i<i,j<m) with m > n. Let 
the rows and columns of A be indexed hy {ii, . . . in} and {ji, . . ■ jn}, respectively. The weight of 
the term corresponding to a in 



det A = y^ sign(cr^ 
a-eSn 



r=l 



2(V+J<.(.)) 



(q), 



(2.25) 



is 2((ii H in) + {ja(i) H Ja{n))) = 2((n H in) + (ji H jn)), which is a constant. 

Next, for an even more general argument, suppose that each term in the n x n determinant 



detA = det{Ep^^{q))i<ij<n= ^ sign{a)Y[Ep^_^^^^{q), 



(2.26) 



a-€S„ 



r=l 



has the same weight, where the n^ subscripts pij are now arbitrary reals. In this setting we take 
the weight of the term in (2.26) corresponding to a to be 



Pl,cr{l) + P2,cr{2) H 1- Pn,cr{n) ■ 



(2.27) 



We then claim that 



Pij -Pij-i=pij -pij-i, for i = 1,2,--- ,n and j = 2,3, • • • ,n. 



and 



Pi,j -Pi-i,j = Pi 



Pi-i, 



for i = 2, 3, • • • , n and j = 1,2,- 



, n. 



(2.28) 



(2.29) 
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To obtain (2.28) and (2.29) first consider the 2x2 submatrix 



(2.30) 



where i, j = 2, 3, • • • , n. Keeping in mind (2.30), there are at least two terms in (2.26) of the form 
B ■ Ep^-^ .J {q)Ep. . (q) and —B • Ep^ ._^ {q)Ep._-^ . {q). Equating the weights of these two terms and 
simphfying, gives 

Pi_ij_i +Pij =Pij_i +Pi_ij, for i,j = 2,3,--- ,n. (2.31) 

By rewriting (2.31) in two ways, we have 

Pi,j -Pij_i =Pi_ij -Pi_ij_i, for i,j = 2,3,--- ,n, (2.32) 

and 

Pi^j-Pi-i^j =Pi^j-i-Pi-i^j-i, for i,j = 2,3,--- ,n. (2.33) 

Equation (2.28) is immediate from (2.32) by fixing j and varying i, while (2.29) is immediate from 
(2.33) by fixing i and varying j. 

By permuting the columns and then the rows of the matrix A in (2.26), and factoring out a 
— 1 if necessary, we can assume that the differences in the right-hand-sides of (2.28) and (2.29) are 
strictly positive. 

It is now clear that if we take the n^ subscripts pi^j to be even integers greater than 2, then 
detvl is ±1 times some n x n minor of a Hankel determinant det(i?2(i+j)(9))i<i,j<m of Eisenstein 
series, with m > n. D 

The simplest application of Proposition 2.4 involves applying the modular forms approach of 
[31, pp. 88-93], as outlined in [5], to express the n x n Hankel determinant Hn ({£'2(;y+i) (*?)}) as 
the product of A"~^(r), a homogeneous polynomial in E^ and Eq, and if needed, E^. One of our 
original motivations for studying the determinants Hn ({-E'2(;y+i) ('?)}) was to note that 

E^ = H[^\{E,^,+,){q)}) = det\E^\, (2.34) 

and recall the n = 2,3,4 cases in (1.6), (1.7), and (2.19), respectively. Since E4 only appears as 
a factor in the left-hand-sides of (2.34) and (2.19), it was natural to split n up into the classes 
(mod 3) given by n = 3r + 1, 3r -|- 2, and 3r -|- 3, with r = 0, 1, 2, • • • . 

Mathematica [37] computations of Hn \{E2(u+i){q)}) up to n = 10, analogous to those of (1.7) 
and (2.19), first led to the discovery of the general form of the following three infinite families of 
formulas in (2.35)-(2.37) below. We prove the evaluations in (2.35)-(2.37) by appealing to the 
standard structure theory of the ring of modular forms in [31, pp. 88-93]. The proof here is a 
detailed rewriting of the original proof supplied by Borcherds in [5]. We have the following theorem. 

Theorem 2.5. Let q := exp(27rir), where r is in the upper half-plane TC. Let A(r) and i?2m = 
E2m{q) be determined by Definitions 1.1 and 1.2, respectively. Then, for \q\ < 1, 

E4A^^{T)-Psrir-i)/2{ElEl)= drHill,{{E2^,+,){q)}), for r = 0,l,2,---, (2.35) 



A'-+\T)-Qri3r-i)/2{ElEl)= erHill^{{E2^,+,){q)}), for r = 0,l,2,..., (2.36) 
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A^'^+\T)-Rr^Sr+i)/2{ElEl)= frHill,{{E2^,+i){q)}), foT r = 0,l,2,..., (2.37) 

where dr, Cr, and fr are constants depending on r, and Pn{x,y), Qn{x,y), and Rn{x,y) are 
homogeneous polynomials in x and y of total degree n (as given above), with integer coefficients, 
whose monomials are those in {x — y)". 

Proof. Let Hn ({£'2(1^+1) (f?)}) denote any of the Hankel determinants in (2.35)-(2.37). Keeping 
in mind that E2r{q) has weight 2r, it is immediate from Proposition 2.4 that each term in the 
Hankel determinant Hn ({£-2(1^+1) ((?)}) is an entire modular form of fixed weight 2n(n +1), as is 
the entire determinant. 

Each Eisenstein series E2r{q) in (1.3) is written in [31, Eqn. (34), pp. 92] as a Maclaurin series 
in q starting with the terms 1 + aiq, with oi 7^ 0. Thus, subtracting the first row from each of the 
others, factoring q out of each of the resulting n — 1 lower rows, and recalling q := exp(27rir), we 

find that Hn \{E2{u+i){q)}) vanishes to order n — 1 at the cusp r = zoo. 

The Maclaurin series in q for A{q) in (1.1) starts with the term q. The function A(g) = A(t) 
is also a cusp form of weight 12 which vanishes at r = ioo. It follows from [31, Theorem 4 (iii), 
pp. 88] that A'^~^{t) divides Hn ({i^2(!/+i) (?)}); and that the quotient is a holomorphic modular 
form of weight 2n(n + 1) — 12(n — 1) = 2(n — 2)(n — 3). 

By Corollary 2 of [31, pp. 89] the quotient Hn ({-B2(;/+i)(9)})/A"~^(r) is a linear combination 
of monomials 

E2E^, (2.38) 

with fixed weight 

4a + 6/3 = 2(n-2)(n-3), (2.39) 

where a and /3 are nonnegative integers. To obtain the left-hand-sides of (2.35)-(2.37) we next 
utilize the weight 2(n — 2)(n — 3) mod 12 to simplify (2.38). 

First, as in [31, pp. 90], let 4a-|-6/3 = 12?tt,, with m a nonnegative integer. This gives 2a + 3(3 = 
6m, which when divided by 2 and 3, respectively, implies that /3/2 = b and a/3 = a are nonnegative 
integers. This gives 

E2E^ = EfEl\ (2.40) 

which in turn leads to (2.36) and (2.37) where n = 3r + 2 and 3r-|-3, respectively. The total degree 
of the polynomials Q and R in (2.36) and (2.37) is immediate from (2.39) by computing 

(4a + 6/?)/12 = a + b={n-2){n- 3)/6, (2.41) 

for n = 3r + 2 and 3r -|- 3, respectively. 

Next, to obtain (2.35), let Aa + Q(3 = 12m + 4. This gives 2a + 3/? = 6m -|- 2. We must have 
a > 1 or else (3 is not an integer. Dividing by 2 implies that (3/2 = 6 is a nonnegative integer. 
Setting (5 = 2b and solving for a gives a = 3(m — 6) -|- 1 := 3c -|- 1, with c a nonnegative integer. 
We now have 

E2E^^ = Ei ■ EfEl\ (2.42) 

which in turn leads to (2.35) when n = 3r + 1. In this case we have 

(4a + Qf3- 4)/12 = c -f 6 = [(n - 2)(n - 3) - 2]/6 (2.43) 

which equals 3r(r — l)/2. D 

It is an interesting open problem to find a concise combinatorial and/or analytical description 
of the coefficients in the polynomials Pir{r-i)/2-, Qr(3r-i)/2, and Rr{3r+i)/2- 
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The r = cases of (2.35), (2.36), and (2.37) are (2.34), (1.6), and (1.7), respectively. The first 
nontrivial case of (2.35) is (2.19). The degree of P3,r{r-i)/2 ™ (2.35) is 3(2), while the degrees of 
Qr(3r-i)/2 and -Rr(3r+i)/2 i^ (2.36) and (2.37), respectively, are the pentagonal numbers r(3r=Fl)/2 
in [32, sequence M1336]. See also [6, pp. 124]. 

Given (2.14), (2.15), (2.20)-(2.22), (2.24), and the modular forms proof of Theorem 2.5, we can 

obtain more complicated infinite families of identities analogous to (2.35)-(2.37) for Xn ({-E'2(!y+i)('?)})- 

Borcherds also observed in [5] that the method of proof of Theorem 2.5 also establishes the 
determinental identities in Theorem 2.3, up to some constant. The space of modular forms of 
the appropriate weight happens to have dimension 1, and is thus spanned by an Eisenstein series. 
The argument through equation (2.39) is the same, with the right-hand-side of (2.39) replaced 
by the sum of the subscripts of the Eisenstein series in the diagonal entries of the given n x n 
matrix, minus 12(n— 1). Call this expression Wi. We then look at Wi mod 12, as before. We find 
that any holomorphic modular form of weight 0,2,4,6,8, or 10 mod 12 is equal to 1, E14, £^4, Eq, 
Es, or Eio times a homogeneous polynomial in E^ and Eg. In the case of Theorem 2.3 we have 
Wi = 4,6,8,10,14. These weights correspond to the right spaces of dimension 1 in our list mod 12, 
the homogeneous polynomial is a constant, and we are done. We have to be careful in the case of 
Wi = 14. Here, we start with Wi = 4a + 6(3 = 12m + 2, for m > 1. We end up factoring out 
E^Eq, with the homogeneous polynomial having total degree m — 1. The other cases are simpler, 
use m > 0, and the homogeneous polynomials all have total degree m. 

The above modular forms proofs of Theorems 2.3 and 2.5, combined with the characterization 
in Proposition 2.4, lead to only a small number of determinental identities in which the space of 
modular forms of the appropriate weight Wi has dimension 1. In particular, as soon as Wi > 16, 
the dimension is > 2, and the homogeneous polynomial in E^ and Eq is no longer a constant. 
All such identities analogous to those in Theorems 1.3, 1.4, and 2.3 are determined by requiring 
Wi < 14. It is not hard to see that this is not possible if n > 6. The remaining finite number 
of possible cases for n < 5 can be checked directly. The identities in Theorems 1.3, 1.4, and 2.3 
cover the basic types that are possible. Moreover, the identities in (1.6) and (1.7) are the only ones 
involving just a power of A, and the identity in (2.24) is unique up to transposition symmetry of the 
5x5 determinant. There are no other analogous 5x5 determinental identities whose left-hand-side 
is of the form £^2r A'^(r). Thus, (2.24) should be very interesting. 

Keeping in mind [35, Eqn. (52.6), pp. 201] it is natural to consider the Hankel determinants 

Hn {{E2{u+i){q)}) in which entries E2[iy-^-i){q) are replaced by unless 2{u + 1) satisfies any of a 
fixed set of congruence conditions. (The "unless" can also be replaced by "whenever"). That is, 
when all entries in certain of the counter diagonals are 0. For example, the condition i?2(!^+i) (q) *~^ 
unless 2{i^ + 1) = (mod 6), leads to interesting determinant evaluations. Similarly, the condition 
i?2(;y+i)(9) ^^ whenever 2(z^ -|- 1) = (mod 4) or 2(z^ -|- 1) = (mod 6) leads to reasonable 
determinants. This is just a small sample of many such possibilities. 

3. The Jacobi elliptic function ns^ and Eisenstein series 

In this section we follow Jacobi's analysis in [16, Section 42] and utilize the Fourier series for the 
Jacobi elliptic function ns^ to write down a formula for the Eisenstein series E2n, for n > 2. We 
then apply [16, Eqn. (2.), Section 36] to put together a simple verification proof of the classical 
formula for A in (1.5). 

We require the Jacobi elliptic function parameter 



z:=2Fi 2^2 J.2 = 2K{k)/n = 2K/Tr, (3.1 

with 

f^ df 
K{k) = K:= = 5 2i^i 



" 1 1 

2 ' 2 
1 


k^ 



(3.2) 
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the complete elliptic integral of the first kind in [21, Eqn. (3.1.3), pp. 51], and k the modulus. We 
also need the complete elliptic integral of the second kind 



E{k) = E:-- 



1 - kH'^ 
l-t2 



dt = ^ 2F1 



1 _ 1 

2 ' 2 

1 



(3.3) 



Finally, we take 



q ■= exp{-TrK{Vl - k^)/K{k)) 



(3.4) 



The classical Fourier expansion for ns^, which first appeared in [16, Eqn. (2.), Section 42; Eqn. 
IV., Section 44], is now given by 



ns^ {u,k) = 1 — 



E 
K 



■ CSC 



2 u 



°° 2t7 

./ ^^COS^. 
y2 / ^ 1 _ r,2n z 



n=l 



(3.5) 



More recent references for (3.5) include [12, Eqn. 27, pp. 913] and [36, Ex. 57, pp. 535]. 

The Fourier expansion in (3.5) may be written as a double sum by first expanding the cos ^^^ 
as a Maclaurin series, interchanging summation, and then simplifying. Next, the Laurent series 
expansion for csc^ - is immediate from differentiating the Laurent series for cot 6 in [7, Ex. 36, 
pp. 88]. This analysis yields 



ns^(ti, k) 



1 



E 



1 ^m^-lo2r?^ — 1 ^.2m — 2 






m=l 


m- z^™ 


(2m -2)! 


t 1 \m — l'y2m—l 


r - ^2^^-1^2.1 

.r=l 


^2m-2 


^2m+2 


(2m -2)! 



(3.6) 



with B2m the Bernoulli numbers defined by (1.4). 
The Laurent series expansion for ns^(u, k) is 



1 



ns' 



!(^,/,)= +^(„S^ 



.{k^ 



u 



2m-2 



m=l 



(2m -2)!' 



(3.7) 



where {ns'^)m{k'^) are polynomials in fc^, with k the modulus. 

In what follows we equate the g's in (3.4) and Definition 1.2. That is 2r = iK{\Jl — k'^)/K{k). 

Keeping in mind Definition 1.2, we find that equating coefficients of n^"*"^, for ?ti > 2, in 
(3.6) and (3.7) leads to a formula for -E2m(5^)- Furthermore, applying the Gaufi transformation 
{q^ V^, k^ j^, K ^ {1 + k)K, z H-> (1 + k)z) from [16, Theorem III, Section 37] to the first 
formula yields a corresponding formula for E2miQ)- We have the following theorem. 

Theorem 3.1. Let z := 2K{k)/iT = 2K/tt, as in (3.1), with k the modulus. Let the Bernoulli 
numbers B2m be defined by (1.4). Take {ns'^)m{k'^) to be the elliptic function polynomials of 
k^ determined by (3.7). Let q be as in (3.4). Take E2m{(l) o-s in Definition 1.2, with 2r = 
iK{Vl~k'^)/K{k). Let m = 2,3,4,- ■ 



E2m{q^) = E2W^t) = l-z^ + 



We then have 

(_l\m-l 



m ■ z 



2m+2 



22 



m—1 



B. 



[ns 



2m 



.{k' 



(3. 



2^2 



E2miq) ^ E2miT) = 1 - (1 + k)' Z 
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+ 22m-i ■ b[^ (ns^Ui^k/il + kf). (3.9) 

The m = 2 and 3 cases of (3.8) are given by 

E4{q'') = z''{l-k'' + k^), (3.10) 

E^iq^) = 2*^(1 + P)(l - 2P)(1 - 1^2). (3.11) 

Equations (3.10) and (3.11) appear in [3, Entry 13(i),(ii), pp. 126]. Equation (3.10) is also 
recorded in [4, Eqn. (12.21), pp. 49], and just under equation (6.) of [16, Section 42]. Jacobi's 
derivation of his 8 squares formula in Section 42 of [16] only required him to go as far as i?4(g^). 

Our verification proof of (1.5) is a consequence of (3.10), (3.11), and equation (2.) of [16, Section 
36] written in our notation as 

Ai/4(2r) ^ Ai/4(g2) ^ ^^3^(1 _ ^2)i/2_ (3^^2) 

The fourth power of (3.12) immediately gives 

A(2r) = A{q^) = j^z^\l - k^fk\ (3.13) 

Equation (3.13) also appears in [3, Entry 12(iii), pp. 124] and [4, Theorem 8.3(iii), pp. 29]. 

Substituting (3.10) and (3.11) into the right-hand-side of (1.5), with q replaced by g^, simplifying, 
and obtaining the right-hand-side of (3.13) now completes the proof of (1.5). The z^^ factored out 
quickly and reduced the proof to a computation involving polynomials of low degree in k"^ . 

Our simple verification proof of (1.5) does not seem to have been written down in the literature 
before. 

Acknowledgment: We would like to thank R. Borcherds for pointing out in [5] how the standard 
structure theory of the ring of modular forms yields the above proof of Theorem 2.5, and also how 
these same techniques may be used to give an alternate proof of Theorem 2.3. 
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